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ABSTRACT: 

The  NASTRAN  finite  element  computer  program  was  used 
to  analyze  the  Mode  1  problems  of  a  thin  strip  fatigue 
specimen  with  a  cracked  central  hole  in  tension,  and  a 
single  edge  cracked  specimen  in  tension  and  bending.  A 
regular  eight  node  isoparametric  quadrilateral  element  was 
degenerated  into  a  quarter  point  triangular  element  which 
approximates  the  stress  singularity  at  the  crack  tip. 
Accurate  values  of  the  Mode  I  stress  intensity  factors 
were  computed  using  this  element. 

The  normalized  stress  intensity  factors, 

Kl/ff0(,ra)1/2*  for  cracks  of  various  lengths  were  computed 
for  the  thin  strip  fatigue  specimen  with  a  cracked  central 
hole  in  tension.  The  normalized  stress  intensity  factor 
increased  near  the  hole. 

A  boundary  condition  comparison  between  a  point  lead, 
a  uniformly  distributed  load,  and  a  fixed  grip- 
displacement  controlled  condition  was  made  on  the 
normalized  stress  intensity  factors  for  the  single  edge 
cracked  specimen  in  tension.  The  stress  intensity  factors 
for  the  point  load  condition  and  the  uniformly  distributed 
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load  condition  were  equal  except  for  very  short  specimen 
lengths.  The  stress  intensity  factors  for  the  fixed  grip 
condition  were  significantly  different  from  the  point  load 
and  uniformly  distributed  load  conditions,  except  for  very 
long  specimens. 

The  normalized  stress  intensity  factors  for  the 
single  edge  cracked  test  specimen  subjected  to  a  pure 
bending  moment ,  four  point  bending ,  and  three  point 
bending  were  also  compared.  The  normalized  stress 
intensity  factors  for  the  four  point  bend  specimen  and  the 
pure  bend  specimen  were  equal  except  for  very  short 
specimen  lengths,  however  the  three  point  bend  specimen 
had  slightly  different  results. 
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1.  INTRODUCTION 

The  finite  element  method  is  useful  for  calculating 
stress  intensity  factors  for  complex  geometries  where 
handbook  solutions  do  not  exist.  This  is  especially 
important  in  alroraft  applications  where  the  geometry  and 
loading  are  usually  complicated.  Conventional  finite 
elements  cannot  adequately  model  the  stress  singularity  at 
the  crack  tipi l]1.  Therefore,  special  methods  must  be 
employed  to  calculate  accurate  stress  intensity  values. 
Since  1970  a  number  of  methods  have  been  developed  to 
compute  stress  intensity  factors  with  relatively  low 
computer  coststl-12].  The  most  popular  approaches  for 
computing  stress  intensity  factors  fall  generally  into 
three  categories,  namely  the  direct  methods,  the  indirect 
methods,  and  the  hybrid  or  cracked  element  method. 

In  the  direct  method,  the  stress  intensity  factors 
are  calculated  directly  by  fitting  the  finite  element 
values  of  the  stresses  or  the  displacements  into  the 
asymptotic  elasticity  relations  for  stress  and 
displacement  around  the  crack  tip.  Usually  special 
elements  12-4]  are  used  around  the  crack  tip  which 

1  Numbers  in  brackets  designate  References  at  end  of 
paper . 
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approximate  the  inverse  square  root  stress  singularity. 
The  advantage  to  this  method  is  its  relative  simplicity. 

The  most  popular  orach  tip  finite  elements  are  that 
of  Tracey[2] ,  Barsoum[3],  Eenshell  and  Shaw[4] .  Both 
elements  approximate  the  stress  singularity  at  the  crack 
tip.  Tracey's  element  must  be  implemented  into  the  finite 
element  program,  however  the  element  in  [3]  and  [4]  is 
obtained  by  degenerating  a  standard  isoparametric 
quadrilateral  element  into  a  triangle,  and  moving  the 
midside  nodes  to  the  quarter  position. 

In  the  indirect  method,  the  stress  intensity  factor 
is  computed  through  its  relation  to  other  quantities,  such 
as  energy  release  rate.  The  most  commonly  used  approaches 
are  the  line  integral  method,  energy  release  method,  and 
the  crack  closure  integral  method. 

The  line  integral  (energy)  method  employs  the  J- 
integral  developed  by  Rice (133, 


where  c  is  an  arbitrary  contour  surrounding  the  crack 

A  A 

tip,  T  is  the  stress  vector  acting  on  the  contour,  u  is 
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the  displacement  vector,  and  V  is  the  strain  energy 
density, 

*  “  1 9  ij  d  6  i  y 

The  line  integral  J  is  related  to  the  stress  intensity 
factors  by 


J  -  [Ki2  +  K22] 

E'  (2) 


where  E'  -  E  (plane  stress) 

E'  -  E/(l-  » 2)  (plane  strain) 

Linear  displacement  finite  elements  can  be  used  with 
this  method  resulting  in  lower  computer  costs,  however 
multiple  analyses  with  the  same  mesh  and  different 
boundary  conditions  must  be  performed  to  uncouple  and 
K 2  in  a  mixed  mode  analysis [5]. 

In  the  energy-release  approach  [6,7,8],  the  stress 
intensity  factors  are  computed  from  the  change  in  energy 
resulting  from  successive  crack  tip  positions  along  the  mesh  . 
The  energy  release  rate,  G,  is  related  to  the  change  in 
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the  total  potential  energy,  np  ,  with  respect  to  an 
increase  of  crach  length,  da,  by 


dn  P  AnP 

da  da  .  (3) 


The  total  potential  energy  of  the  finite  element 
model  is 

np  -  (1/2){u}T  [K]  {u}  -  {u}T  {f}  (4) 

where  {u}  is  the  displacement  vector,  [K]  is  the 
stiffness  matrix,  and  { f >  is  the  force  vector.  The  energy 
release  rate  can  be  computed  by  incrementing  the  crack, 
length  by  a  small  amount  da.  then 

drrp  -  (  np  +  dnp)  -  np  -  (i/2){u}T  t  dK]  {u}  cs) 

Extending  the  orach  by  da  affects  the  elements 
around  the  orach  tip  only,  therefore  [  dK]  need  only  to  be 
evaluated  for  those  elements. 
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The  energy  release  rate.  G  ,  is  computed  from 
equation  (3)  and  the  stress  intensity  factors  are  related 
to  G  by 


[*12  +  K2S]  (6) 

G  -  - 

E' 


where  E  -  E  (plane  stress) 

E‘  -  E/(l-  v2)  (plane  strain). 

The  crack  closure  integral  method  permits  both  K]_ 
and  Kg  to  be  computed  from  a  single  analysis  [9).  This 
method  is  based  on  Irwin's  concept  that  if  a  crack  extends 
by  a  small  amount  da,  then  the  energy  absorbed  in  the 
process  is  equal  to  the  work  required  to  close  the  crack 
to  its  original  length.  The  energy  release  rate  equations 
for  mode  I  and  II  loading  are: 

/"da 

Gi  -  lim  1/(2  da)  I  ffyy(  da  -r,0)v(r,  *  )dr, 
d  a-*  o  Jo 
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G2 


-  lim 
da— »o 


1/(2  ia) 


ffxy(  Aa  -r,0)u(r,  *  )dr. 


(?a.b) 


The  integrals  physically  represent  the  amount  of  work 
required  to  close  the  orack  by  a  small  amount  a  a.  In 
terms  of  the  finite  element  representation,  this  work  is 
equal  to  one-half  the  product  of  the  forces  at  nodes  c  and 
d  (Figure  1)  and  the  distance  (vc  -  vd)  required  to  close 
the  nodes.  The  equations  for  G^  and  G2  become 


Gx  -  11m  Fc-(vc-vd) 

aa-MD  2Aa 

(8a, b) 

G2  -  Hm  _1_  Tc-  (Uq-u^) 

Aa-*0  2Aa 

where  Fc  and  Tc  are  the  X  and  Y  forces, 
-respectively,  required  to  hold  nodes  c  and  d  together. 

This  approach  works  best  with  linear  displacement 
elements . 

The  Hybrid  element  approach  utilizes  cracked  elements 
containing  an  exact  mathematical  solution  to  the  stress 
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singularity  within  them  [10,11].  The  stress  intensity 
factors  are  determined  from  nodal  point  displacements 
along  the  periphery  of  the  cracked  element.  These  elements 
are  either  incorporated  into  the  finite  element  mesh  or 
are  solved  separately  using  the  finite  element  results. 
Hybrid  elements  yield  highly  accurate  solutions  with  very 
coarse  grids. 

The  objective  was  to  choose  a  simple  method  for  two- 
dimensional  analysis  that  can  be  used  with  the  NASTRAN 
Finite  Element  Computer  Program,  without  requiring  any 
computer  programing.  For  this  purpose,  the  degenerated 
isoparametric  element  Introduced  by  Barsoum  and  Henshav  et 
al.  is  the  most  attractive  approach.  This  element  used 
with  the  direct  method  can  calculate  both  and  K2  a 
single  analysis. 
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2.  THEORY 

2.1  Linear  Elastic  Fracture  Mechanics 

Consider  a  flat  orach  in  a  two  dimensional  linearly 
elastic  solid  under  ln-plane  loading  shown  in  Figure  2a. 
The  stress  and  displacement  states  at  a  small  distance  P 
from  the  orach  tip  described  by  the  polar  coordinates  T,e 
are  (Figure  2b)  given  by  [14] 


ffxx(r.  e  )  - 


(2  »  r)1/2 
*2 


(2  t  r)1/2 


-cos( 0/2) |^1  -  sia(0/2)sln(30/2) J 
— sin(0/2)[2  +  cos(0/2)cos(30/2)l 


'yyCr,  e  )  - 


*1  r  i 

cos(  6  /2)  1  +  sln(  9  /2)sln(3  6  /2) 


“(2  Vrj'^a 
*2 

(2  T r)1/2 


Sin(  0/2)cos(  6  /2)cos(3 6  /2) 


Si 

*xy(r,  (9  )  -  _ _  sln(  0/2)cos(  0/2)cos(3  0  /2) 

(2  t  r)i/a 
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+ _ _  cosC  0/2)  1  -  sin(  0/2)sin(3  d  /2) 

(2  *  r)  i/a  L 

C9a-d) 

tfzz  "  0  (plane  stress) 

*zz  “  •'C  ®xx  *  *yy)  (plane  strain) 


u(r. 

,  Kl 
e  )  -  _ 

2G 

.  (r/2  «■  )1/2  cos(  0 / 2 ) [^  k-  1  +  2sln2(0/2)J 

*2 

+ 

.  (r/2*  )1/2  sin(0/2)  Fk+  1  +  2COS2(0/2)1 

2G 

L  J 

v(r. 

,  *1 

•  )  -  _ 

2G 

(r/2*)1/2  sin(0/2) 

£  K+  1  -  2COS2(0/2)J 

+ 

1  W 

1  to 

(r/2  r )1/2  cos( 0/2) 

1  -K  t-  2sln2(0/2)j 

2G 

(lOa-b) 


where  G  Is  the  shear  modulus,  *  is  Poisson's  ratio 

and 
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k-  (3-4 *0  (plane  strain) 

k  -  (3-  0/(1+  O  (plane  stress) 

The  constants  and  are  the  stress  intensity  factors 
associated  with  the  loading  modes  shown  in  Figure  3. 
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2.2  Isoparametric  Finite  Element 

The  essential  idea  behind  the  8-node  isoparametric 
element  is  that  a  simple  square  element  defined  in  a  local 
coordinate  system  may  be  transformed  into  a  curve-sided 
element  in  the  global  Cartesian  Coordinate  System,  Figure 
4.  For  an  element  to  be  isoparametric  the  interpolation 
functions  defining  the  geometry  and  the  displacement  field 
are  the  same. 

The  interpolation  functions  of  the  8-node  quadratic 
isoparametric  element  (Serendipity  Type)  in  the  local  £  -  n 
coordinate  system  are  [153: 

*1  -  ( 1/4) ( 1  -  £  )(1  -  n)  -  (i/2)(N5  +  Nq) 

N2  -  ( 1/4) ( 1  +  £  )(1  -V)  -  (1/2)(N5  +  N6) 

N3  -  ( 1/4) (1  +  £  )(1  +0  -  (l/2)(»e  +  N7) 

K4  -  (1/4)(1  -  £  )(1  +  V  )  -  (1/2)(N7  +  Ng) 

N5  -  (1/2X1  -  £  2)(1  -  n  ) 

N6  -  ( 1/2) ( 1  -  «2)(1  +  5  ) 

N7  -  (1/2)(1  -  £  2)(1  + n  ) 

NS  -  (1/2)(1  -  n2)(l  -  £  )  (lla-h) 
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where  the  interpolation  functions  have  the  property, 
Ni  -  1  at  node  i, 

Ni  -  0  at  every  other  node. 

The  napping  which  transforms  the  square  element  from 
local  coordinates  to  global  X-Y  coordinates  is: 


X 


8 


£ 

1-1 


N±(  S  .  )x± 


Y 


8 

E  »i(  (  ,  i  )y± 

i-i 


(12a, b) 


where  1  corresponds  to  the  node  i  whose  coordinates 
are  (X^Y^)  in  the  X-Y  system  and  (  £j_,  *  in  the  £  -  v 
system. 

The  displacements  are  represented  through  the 
interpolating  functions  by: 
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U(  S  ,  1  ) 


8 

S  i  ,  1  )Ui 

i-1 


CNiHui) 


v(  $  ,  V  ) 


N±( 


*7  -  [Nj_]{v^} 


lu 

-  [Hi] 

U1 

(v 

vi 

(13a, b.o) 


where  uj.  and  are  the  calculated  displacements  at 
node  1. 

The  strains  can  be  calculated  from  the  elasticity 
strain-displacement  relationship  and  equations  (13a-c). 
The  strain-displacement  relationship  is 


(<zx 

-  HI 

d/dx 

0 

( *77 

CD]  - 

0 

*/*7 

'2*xy 

_d/dy 

d/dx 

(14) 
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Substituting  equation  (13c)  into  equation  (14),  the 
strains  can  be  represented  through  the  interpolation 
functions  and  the  nodal  displacements  by 


(  *  }  -  tDNNi] 

ui 

Ui 

-  [B3 

vi 

Vi 

[B] 


dNi/dx  0 

0  dtfj/dy 
dNj/dy  dNi/dx 


(15) 


Equation  (15)  requires  differentiation  of  the 
interpolation  functions  with  respect  to  X  and  Y.  Since  the 
interpolation  functions  N^(  S  ,  'i  )  are  functions  of  £  and  n 
the  differentiation  may  be  accomplished  through  the  chain 
rule , 
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dNi/a*  j 

ax/as 

ay/a* 

aN*/  ax  j 

aNj/  a  x  j 

}  m 

-  [j] 

I  , 

aNi/an) 

dx/a  n 

m 

ay/an 

aNi/ay)  I 

aN±/  ay) 

(16) 

where  [J]  Is  the  Jacobian  matrix. 

Thus , 


( aN±/ax 

( a^/ai 

-  ur1 

( a^/ay 

1  ( aNi/av 

(IT) 


Equation  (17)  can  be  used  to  evaluate  the  terms  in 
equation  (15). 

For  a  homogeneous,  isotropic  linear  elastic  body  in 
plane  stress,  the  stresses  are  calculated  from  the  strains 
using  Boohe's  Law. 
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ff-rr 

i - 

o 

i — i 

i _ 

{ *  }  - 

*yy 

-  [C]  {« } .  [Cl  -  E/(l-  *  2) 

"  1  0 

1 

axy 

0  0  (1/2) ( 1-  »  ) 

(18) 


t  E  -  Young ' s  Modulus 

►  v  m  Poisson's  Ratio 
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2.3  Quarter-Point  Triangular  Crack  Tip  Element [3, 4] 

Figure  5  shows  a  triangular  element  generated  by 
collapsing  side  1-4  of  the  quadrilateral  element  in  Figure 
4  and  moving  the  mid-side  nodes  adjacent  to  the  crack  to 
the  quarter  position.  Side  1-4  was  arbitrarily  chosen  to 
be  the  collapsed  side.  The  crack  runs  along  the  negative 
X-axis . 

Along  side  1-2  of  the  triangle,  n  -  -1  and  the 
interpolation  functions  are: 

N3  -  N4  -  N6  -  N7  -  NS  -  0 

NX  -  (-1/2)  5(1  -  5  ) 

N2  -  (1/2)  5(1  +  5  ) 

»5  -  Cl  -  « 2)  (21) 

From  equation  (12) 

X  -  (-1/2)5(1  -  5  )*!  +  (1/2)  5(1  +  (.  )x2  +  (1  -  52)x5 

(22) 
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Choosing  X^  »  0 ,  X2  m  ^5  “  L/4,  then 

X  -  (l/2)£ (1  +  £  )L  +  Cl  -  £2)(L/4) 

(23) 

Thus, 

£-  (-1  +  2(x/L)1/2  )  (24) 

From  equation  (13),  the  displacement  along  side  1-2  is, 
u  -  (-l/2)£(l-  £ )u!  +  (l/2)£(l+  £  )u2  +  (l-  £2)u5 

(25) 

or  in  terms  of  X  is, 

u  -  (-1/2) (-1  +  2(x/L)1/2)[2  -  2(x/L)1/2]u1 

+  ( 1/2) (-1  +  2(x/L)1/2)[2(x/L)1/2]u2 

+  (  4(x/L)1/2  -  4(x/L)  )U5 

(26) 

The  displacement  relation  in  equation  (26) 
approximates  the  r1/2  behavior  near  the  cracb  tip  as 
equation  (10a). 

The  strains  in  the  X-direction  are, 
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*xx  -  du/ax  -  (-1/2)  [3/(xL)1/2  -  4/lJux 

+  (1/2) [  -l/(xL)1/2  +  4/Lju2  +  ^2/ (xL)1/2  -  4/lJu5 

(2?) 

Equation  (27)  has  the  inverse  square  root  singularity 
at  the  orach  tip  as  equation  (9).  Therefore,  the  Linear 
Elastic  Fracture  Mechanics  Equations  are  approximated  in 
the  quarter-point  orach  tip  element. 
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3.  NUMERICAL  ANALYSIS 

The  NASTRAN  Finite  Element  Computer  Program  was  used 
to  analyze  the  following  problems:  (1)  a  central  cracked 
plate  test  oase  used  to  study  the  accuracy  of  the  orach 
tip  elements,  (2)  a  thin  strip  fatigue  specimen  with  a 
cracked  central  hole,  (3)  a  boundary  condition  comparison 
on  a  6ingle  edge  cracked  specimen,  and  (4)  a  comparison  of 
single  edge  cracked  bending  specimens.  In  all  problems, 
the  following  assumptions  are  applied: 

1.  the  specimens  are  made  of  7075-T6  Aluminum  which  is 
homogeneous, isotropic,  and  linear  elastic, 

2.  the  specimens  are  thin  plates  under  in-plane  Mode  I 
loading  conditions, 

3.  the  problems  are  in  the  plane  stress  condition. 

Collapsed  quarter  point  quadrilateral  elements  were  used 
around  the  crack  tip  while  the  8-node  isoparametric 
quadrilateral  elements  were  used  in  the  rest  of  the  mesh. 
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3.1  The  Centrally  Cracked  Plate  Test  Case 

The  crack  tip  elements  were  Initially  tested  on  a 
uniformly  loaded  center  cracked  plate,  shown  in  Figure  6, 
to  study  their  accuracy.  Only  one  quarter  of  the  problem 
was  analyzed  because  of  symmetry.  The  following  boundary 
conditions  apply  (Figure  6a) : 

1.  cr^Cb.y)  -  ffxy(b.y)  -  0. 

2.  CTjjy(X,h)  ”  0,  tT  yy(X,h)  *  '  Q 

3.  ffxy(0 ,y)  -  0.  u(0,y)  -  0, 

4.  ffxyCx.O)  -  0,  0  <  x  <  b, 

ffyy(X.O)  -  0.  0  <  x  <  a. 

v(x,0)  -  0,  a  <  x  <  b. 

A  total  of  66  elements  and  229  grid  points  were  used 
in  the  mesh  shown  in  Figure  7.  The  quarter-point  element 
side  length  to  half  crack  length  ratio,  L/a  was  L/a  w.io 
The  normalized  Mode  I  stress  Intensity  factor  Hi/  c  i**).1 
calculated  from  the  finite  element  results,  was  1.146 
which  is  l.l  percent  different  from  the  theoretical  value 
of  1.159  obtained  from  [16), 

Ki/  a0(r  a)1/2  -  [  sec(*-a/2b)  ]1/2 

(28) 
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An  accurate  value  of  the  Mode  I  stress  intensity  factor 
was  computed  with  a  relatively  coarse  mesh. 

The  following  guidelines  should  be  applied  when  using 
the  collapsed  quarter  point  quadrilateral  elements  for 
fracture  mechanics  analyses: 

1.  the  grid  points  at  the  crack  tip  (on  the  collapsed 
side)  must  be  rigidly  constrained, 

2.  the  sides  of  these  crack  tip  elements  must  be  kept 
straight  to  avoid  error  [l?j, 

3.  the  L/a  ratio  should  be  kept  small,  L/a  i  .2,  based 
on  this  test  case  and  on  the  results  from  [18], 

4.  the  finite  element  displacement  results  should  be 
used  to  calculate  the  stress  intensity  factors 
(direct  method)  since  they  are  more  accurate  than 
the  stresses, 

5.  the  crack  opening  displacements  v(r,  «■  )  should  be 
used  to  compute  and  u(r,  ~  )  should  be  used  for 
K2,  (equations  (10a,b)). 
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3.2  Thin  Strip  Fatigue  Specimen  With  A  Cracked 
Central  Hole 

The  thin  strip  fatigue  specimen,  Figure  8,  was  used 
to  study  a  hole-crack  interaction  problem.  The  particular 
dimensions  in  Figure  8  represent  the  lower  front  spar 
region  of  the  A-4F  Blue  Angel  aircraft  where  fatigue 
cracks  have  occurred.  The  specimen  was  gripped  at  the 
short  ends  and  displaced  in  the  Y-direction.  The  problem 
was  analyzed  in  two  parts:  (l)  fatigue  specimen  without  a 
crack  and  (2)  fatigue  specimen  with  a  cracked  central 
hole . 
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3.2.1.  Fatigue  Specimen  Without  A  Crack 

A  static  stress  analysis  was  performed  on  the 
specimen  using  finite  elements  to  determine  regions  of 
high  stress  and  the  elastic  stress  concentration  factor, 
Kt.  Due  to  symmetry,  only  one  quarter  of  the  problem  was 
considered.  The  following  boundary  conditions,  shown  in 
Figure  9,  apply: 

1-  -oxy(b-y)  -  o. 

2.  -xy(0»y)  “  0, 

u(0,y)  -  0, 

3.  xyCx.O)  -  0. 

v(x , 0)  -  0, 

and  at  the  fixed  grip  end, 

4.  u(x,d)  -  0,  v(x , d)  -  v0. 

The  nominal  stress  along  the  fixed  grip  edge,  ~0,  was 
calculated  by 

?1 

_  (29) 

A 
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where  are  the  nodal  reaction  forces  Cat  y-d)  in  the  Y- 
direction  from  the  fixed  grip,  and  A  is  the  cross- 
sectional  area  at  the  fixed  grip  end. 

The  finite  element  mesh.  Figure  10,  consisted  of  162 
elements  and  557  grid  points.  The  maximum  stress  occurred 
at  the  corner  of  the  large  central  hole  on  the  X-axis  and 
the  elastic  stress  concentration  factor  based  on  gross 
area,  K^g.  is  3.99.  The  tangential  stress  distribution, c e e 
is  provided  in  Table  1,  Figure  11  for  the  central  hole  and 
in  Table  2,  Figure  12  for  the  small  hole. 
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3.2.2  Fatigue  Specimen  With  A  Cracked  Central  Hole 

The  effect  of  a  crack  of  various  lengths  in  the  large 
central  hole  of  Figure  8  was  considered.  The  crack  is 
located  along  the  positive  and  negative  X-axis.  The 
boundary  conditions  for  this  problem.  Figure  13,  are: 


1. 

axxO>,y)  m  ~ xy 

(b.y)  -  0, 

2. 

o 

i 

>> 

o 

ft 

D 

u(O.y)  -  0, 

3. 

Q 

3 

s-\ 

H 

O 

w 

1 

O 

0  <  x  <  b. 

o 

1 

o 

H 

t> 

0  <  x  <  a, 

v(x, 0)=0. 

a  <  x  <  b. 

and  at  the  fixed  grip  end, 

6.  u(x,d)  -  0 
v(x,d)  -  v0. 

The  nominal  stress  along  the  fixed  grip  edge  was 
calculated  from  equation  (29). 

Figures  14a, b,c  shov  a  typical  mesh  used  in  the 
cracked  specimen  analysis.  A  refined  mesh  was  used  in  the 
region  containing  the  crack.  The  number  of  elements  ranged 
from  228  to  287,  and  the  number  of  grid  points  ranged  from 
765  to  956,  more  elements  were  used  as  the  crack  length 
Increased. 
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Figure  14c  shows  an  enlarged  view  of  the  orach  tip 
region.  The  L/a  ratio  vas  kept  smaller  than  L/a  i  .15. 

This  crack  tip  region  remained  the  same  and  was  moved 
along  the  positive  X-axis  as  the  crack  length  increased. 

The  normalized  stress  intensity  factors  for  various 
lengths  are  provided  in  Figures  15,16  and  Table  3.  A  crack 
length  of  .375  inches  corresponds  to  (c+a)-b  when  the 
crack  reaches  the  side,  hence  the  normalized  stress 
intensity  factors  Increase  toward  infinity.  Even  though 
the  normalized  stress  intensity  factor  decreases  before 
starting  to  increase  with  the  increasing  crack  length 
(Figure  15),  the  absolute  value  of  the  stress  Intensity 
factor  is,  as  expected,  a  monotonically  Increasing 
function  of  the  crack  length  (Figure  16). 

The  limit  of  the  normalized  stress  Intensity  factor 
as  the  crack  length  approaches  zero  in  Figure  15  should 
approach  the  theoretical  solution  to  the  semi-infinite 
plate  with  an  edge  crack.  As  the  crack  length  becomes 
infinitesimal,  the  radius  of  curvature  of  the  hole,  and 
the  length  and  width  of  the  plate  appear  to  be  infinitely 
large.  The  stresses  around  the  infinitesimal  crack  are  no 
longer  affected  by  the  finite  geometry.  The  stress 
intensity  faotors  for  a  uniformly  loaded  semi-infinite 
plate  with  an  crack  is  (16]: 
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Li  -  1.1215*  maz(  it  a) 1/2 .  Z2  -  0. 


(30) 


where  0  Bax  is  the  stress  applied  normal  to  the 
crack . 

On  the  fatigue  test  specimen,  the  stress  in  the 
corner  of  the  hole  where  the  crack  is  located  is : 

c  mam  “  3.99  (31) 

Substituting  equation  (31)  into  equation  (30),  the  stress 
Intensity  factor  becomes, 

*1  -  4.4748  ~0(  ra)1/2.  (32) 

Normalizing  ohe  stress  intensity  factor  in  equation  (32), 
&1 I  -  o  (  »  a) 1/2  -  4.4748,  (33) 

hence  the  curve  in  Figure  15  should  approach  the 
value  of  4.4748  near  a  -  0.  As  seen  from  the  dashed  line, 
the  curve  is  accurate. 
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The  normalized  stress  intensity  factor,  (Kj/  ^0) 
curve  in  Figure  16  should  go  to  zero  as  the  crack  length 
approaches  zero.  The  stress  Intensity  factor  in  equation 
(30)  can  he  written  in  the  form, 

Kj/  °o  “  constant •  (a)1/2  (34) 

Equation  (34)  has  a  parabolic  shape  and  as  the  crack 
length  approaches  zero,  Kj,/  °0  goes  to  zero.  The  finite 
element  results  do  agree  with  equation  (34)  near  a  -  0, 
and  therefore  the  normalized  stress  intensity  factor  curve 
should  follow  the  dashed  line  in  Figure  16. 
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3.3  An  Applied  Load  Boundary  Condition  Comparison  On 
A  Single  Edge  Cracked  Test  Specimen 

Tlie  object  of  this  problem  was  to  determine  the 
effect  of  different  loading  conditions  on  the  stress 
intensity  factor  for  a  single  edge  cracked  test  specimen 
shown  in  Figure  17.  Three  types  of  tensile  loading  were 
considered:  (1)  a  point  load  assumed  to  be  equivalent  to  a 
pin  loaded  experiment,  (2)  a  uniformly  distributed 
load, (3)  an  external  load  applied  through  a  rigid  or  fixed 
grip.  The  length  of  the  specimen  was  varied  to  determine 
at  what  length  (invoking  the  St.  Venant ' s  principle) 
statically  equivalent  loads  applied  in  different  ways 
would  produce  the  same  Xj  values. 

Only  half  of  the  problem  was  analyzed  due  to 
symmetry.  The  boundary  conditions  for  the  specimen  with  a 
point  load,  P,  (Figure  18a)  are: 

1-  ^xx(w-y)  "  'xy(v.y)  -  0, 

2 .  3  xx(O.y)  “  Jxy(^>y)  “  o, 

3.  -  xy(x,0)  -  0,  0  <  x  <  w, 
cyy(X,0)  “0,  0  <  x  <  a, 

v(x,0)  -  0,  a  <  x  <  w, 

4.  0  yy(x,h)  -  P6(x  -  w/2)  ,  ~  xy(x,h)  -  0. 
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and  the  nominal  stress  at  the  loaded  end  (y-h) .  c  o>  is 
given  by 

0  o  -  P/A 

(35) 

where  A  is  the  cross  sectional  area. 


Boundary  conditions  1-3  apply  for  the  specimen  with 
a  uniformly  distributed  load  (Figure  18b),  along  with  the 
condition. 


4.  “  yy(x,h)  “  ^  Q  ,  C  xy(x>k)  “  0. 

In  the  finite  element  analysis,  :  0  was  calculated  from. 


~  o  “  - 

A 


(36) 


where  Pi  are  the  node  point  forces  in  the  Y -direction. 

For  the  fixed  grip-displacement  controlled  problem 
(Figure  18c),  boundary  conditions  1-3  apply,  along  with 
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the  condition, 

4.  u(x,h)  -  0,  v(x,h)  -  Vo,  0  <  x  <  w, 

and  o  o  is  calculated  from  equation  (29). 

Typical  finite  element  meshes  for  this  analysis  are 
shown  in  Figures  19a,b,c.  The  crack  length,  a,  and 
specimen  width,  w,  were  constant  with  a/w  -  .5,  while  the 
specimen  length,  h,  was  varied  between  h/w  -  .1  and  h/w  - 
100.  The  ASTKC19]  dimensions  for  this  specimen  are  a/w  - 
. 5  and  h/w  -  1.5. 

The  normalized  stress  intensity  factor,  Ki/-0(r-a)1/2 , 
is  shown  for  different  half-length  to  width  ratios  (h/w) 
in  Figures  20a, b  and  Table  4.  The  normalized  stress 
intensity  factor  of  2.836  from  Table  4  is  0.2  percent 
different  than  the  theoretical  value  of  2.830  for  the 
uniformly  loaded  strip  with  an  edge  crack  [16].  For  h/w  < 
1.25,  the  point  and  distributed  loading  conditions  have 
different  normalized  stress  Intensity  values,  but  for  h/w 
Z  1.25  they  are  exactly  the  same.  The  fixed  grip  condition 
produces  very  different  results  from  the  point  and 
distributed  load  conditions  for  h/w  <  20.  For  very  small 
h/w  ratios,  the  normalized  stress  Intensity  factors  for 
the  point  and  distributed  loading  conditions  significantly 
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increase,  while  the  normalized  stress  intensity  factor  for 
fixed  grip  condition  continues  to  decrease. 

Deformed  plots  of  the  single  edged  cracked  specimen 
under  the  different  applied  loading  conditions  are 
provided  in  Figures  21a,b,o,d,e,f  for  two  different 
specimen  lengths.  The  deformed  plots  in  Figures  21a, b,o 
show  that  the  specimen  rotates  for  the  point  and 
distributed  load  cases,  however  the  rigid  fixed  grip 
prevents  rotation.  Also,  the  crack  opening  displacements, 
v(r,  r  )  are  larger  for  the  point  and  distributed  load 
conditions,  than  for  the  fixed  grip  conditions  which  is 
why  the  stress  intensity  factor  is  lower  for  the  fixed 
grip  in  Figures  20a, b.  The  crack  opening  displacements 
change  as  the  specimen  length  decreases  which  support  the 
Increase  in  the  normalized  stress  intensity  factors  for 
the  point  and  distributed  loads  and  the  decrease  for  the 
fixed  grip  for  short  specimen  lengths. 

The  normalized  transverse  relative  displacement, 

Au-t/( 30h/E) ,  is  shown  in  Figures  22a, b  and  Table  5.  The 
transverse  displacement,  au-^,  is  defined  as, 

aut  -  (aui)  +  |au2|  (37) 

where  au^  and  au2  are  the  displacements  shown  in  Figure 
21a, b,c.  The  normalized  transverse  relative  displacements 
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for  the  fixed  grip  condition  are  smaller  than  for  the 
distributed  and  point  load  conditions  because  Aug  “  0  for 
the  fixed  grip  conditions. 

The  normal  stress  distribution  in  front  of  the  crack 
tip,  Oyy(x,0),  for  various  h/v  values  is  provided  in 
Figures  23a, b,c  and  Tables  6a,b,c.  For  the  point  and 
distributed  load  conditions,  the  normal  stress 
distribution  in  front  of  the  crack  tip  changes  very  little 
as  h/w  decreases  from  100  to  1.5,  however  the  distribution 
changes  significantly  as  h/w  decreases  from  1.5  to  0.1. 
Significant  changes  in  the  stress  distribution  occur  for 
the  fixed  grip  condition  for  each  h/w  value  in  Table 
6c.  For  specimens  of  long  length  (h/w  -  100),  there  are 
compressive  stresses  near  the  edge  of  the  specimen  [near 
x/w  -  1.0) . 

The  fixed  grip-displacement  controlled  problem  was 
studied  further  because  of  its  complexity.  The  finite 
element  results  were  compared  to  the  two  limiting  analytic 
solutions:  Cl)  the  tearing  of  an  infinitely  wide  strip 
[16]  for  a  very  short  specimen.  Figure  24,  and  (2)  the 
semi-infinite  strip  with  a  fixed  end  [20]  for  a  very  long 
specimen.  Figure  25. 
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The  tearing  of  an  infinitely  wide  strip  problem  is 
shown  in  Figure  24.  An  external  load  is  applied  through 
rigid  clamps  held  at  *yy  - *o  when  the  material  begins  to 
fraoture.  The^e  are  three  distinguishable  zones  in  the 
strip,  zone  I,  an  unloaded  region  where  the  strain  energy 
density  Of)  equals  zero,  zone  II,  where  t  is  unknown,  and 
zone  III,  where  W  -  constant,  ‘yy  - *o  ,  *xx  -  ‘ xy  -  0 . 

The  stress  Intensity  factor  is  given  by  [16] 

Kiv  -  «0  E  (h/{L-  w2))172  (38) 


The  strain  «0  is  given  by, 


(39) 

where  V0  is  the  displacement  at  the  fixed  grip  end  and  E 
is  Young's  Modulus.  Substituting  equation  (39)  into 
equation  (38),  the  stress  intensity  factor  can  be 
expressed  as 
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*1V 


E  v0 

[h(l-  ».  2)  ]  1/2  _ 


(40) 


A  comparison  between  the  normalized  stress  intensity 
factors  for  the  single  edge  cracked  test  specimen  under 
fixed  grip  conditions,  and  the  infinitely  wide  s;rip  with 
a  fixed  end,  is  provided  in  Figures  26a, b  and  Table  7.  The 
stress  intensity  factors  are  normalized  by  the  right  hand 
side  of  equation  (40) ,  which  means  that  the  normalized 
stress  intensity  factor  for  the  infinitely  wide  strip, 

K]_v,  has  a  constant  value  of  one  for  all  h/w  values.  The 
normalized  stress  intensity  values  for  the  single  edge 
cracked  specimen  are  less  than  five  percent  different  from 
the  infinitely  wide  strip  for  h/w  values  less  than  five. 
The  results  for  the  specimen  diverge  from  the  infinitely 
wide  strip  solution  for  large  h/w  values  because  of  the 
long  specimen  length. 

The  normal  and  shear  stress  distribution  along  the 
fixed  end  of  the  single  edge  cracked  specimen  is  shown  in 
Figures  27  and  Table  8  for  various  h/v  values.  For  the 
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short  specimen  (h/v  -  .1)  in  Figure  27a,  the  unloaded 
region  (zone  I)  can  be  seen  near  x/v  -  0. 

The  normal  stress  at  the  fixed  grip  end,  °  yy(x,h), 
(Figure  27a  and  Table  8)  and  the  normal  stress  in  front  of 
the  crack  tip,  a  yy(x,0),  (Figure  23c  and  Table  6c)  are 
not  constant  near  x/v  -  1  (zone  III)  which  means  that  the 
strain  energy  density  is  not  constant  throughout  zone  III. 
This  supports  the  differences  in  the  normalized  stress 
intensity  factors  between  the  specimen  and  the  Infinitely 
wide  strip  in  Figure  26a, b  and  Table  7. 

The  semi-infinite  strip  solution  with  a  fixed  end  by 
Gupta  [20],  shown  in  Figure  25,  was  used  to  study  the 
stresses  around  the  fixed  grip.  The  solution  shows  that 
the  shear  and  normal  stresses  are  singular  at  the  corners 
of  the  fixed  end.  The  stress  intensity  factors  at  the 
corners  are  defined  by, 

51  -  lim  (2)1/2*r“  ~  ^(r ,h) 

r-*Q 

52  -  lim  (2)1/2Ta  o  (r.h)  (41) 

r-*0 


where  r  is  a  variable  starting  at  the  corners  of  the  fixed 
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edge  shown  in  Figure  25.  The  characteristic  equation  to 
determine  a  is  given  by 

2/CCOS7T a  -  K2  +  1  -  4(  <*  -  l)2 


(42) 

where  k-  3-4*  (plane  strain) 

Km  (3  -v  )/(l+  *  )  (plane  stress). 

The  exponent  of  the  stress  singularity,  a  ,  depends  on 
Poisson's  ratio.  The  results  of  [20]  also  show  that  the 
stress  intensity  factor  Ki  depends  on  K2  and  the  ratio 
K2/K1  is  given  by 


Ko  (  k  +  1 )  Sin7T  a 

—  "  -  - 5  (43) 

Kl  f(<-l)(COS7T  a  +  1)  -  2(K+1)(<*-1)  +  4(cr-l)'i] 


For  the  plane  stress  conditions  and  *  -  .33,  <*  -  .2539 

from  equation  (42)  and  K2/Ki  -  -.255  from  equation  (43;. 

The  finite  element  results  for  the  single  edge 
cracked  specimen  (h/w  -  100)  were  compared  to  Gupta’s 
results.  The  stresses  on  the  fixed  edge  of  the  specimen, 
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shown  in  Figure  27e,  compare  well  with  the  results  from 
[20] .  except  at  the  corners  where  the  stress  singularities 
exist.  The  finite  elements  did  not  provide  an  accurate 
representation  of  the  normal  and  shear  stress 
singularities  near  x/w  -  0  and  x/w  -  1,  as  shown  in 
Figures  27a,b,c,d,e  ,  since  the  stress  variation  in  the 
elements  is  linear.  The  actual  stress  distribution  should 
follow  the  dashed  curve.  The  6hear  stress  in  Figure  27e  is 
an  odd  function  in  X,  as  it  was  in  [20]. 

From  equations  (41),  the  stress  intensity  factor 
ratio  can  be  expressed  as. 


Kg 

*1 


r  xy  (r,h) 

lim  - - - 

r~>0  ~  yy  (  r ,  h ) 


(44) 


The  limit  in  equation  (44)  was  approximated  from  the  r  -  0 
value  of  the  finite  element  stress  ratios  shown  in  Figure 
28  and  Table  9.  For  side  A,  Kg/Ki  -  -.2154  and  for  side  3, 
K2/Ki  -  .2107  which  are  15.5  and  17.4  percent  different 
from  the  analytical  value. 
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3.4  A  Comparison  of  Single  Edge  Cracked  Bend 
Specimens 

A  comparison  of  the  following  types  of  single  edge 
cracked  bending  specimens  is  made:  (1)  pure  bending 
applied  through  a  linear  stress  distribution  normal  to  the 
edge  (Figure  29a),  (2)  four  point  bending  (Figure  29b), 
and  three  point  bending  (Figure  29c).  Due  to  symmetry, 
only  one  half  of  the  problem  was  analyzed. 

The  boundary  conditions  for  the  pure  bend  specimen 
are  (Figure  30a): 

l*  r  xxO.y)  -  =xy(v.y)  -  o, 

2.  -  3xy(°-y)  ■  °* 

3.  -^(x.O)  -  0,  0  <  x  <  w, 

~yy(X,0)  “0,  0  <  x  <  a, 

v(x,0)  -  0,  a  <  x  <  w, 

4.  3  yy(x,h)  “  30  ,  '  xy(x.h)  =  0, 

where  '  0  is  a  linear  stress  distribution  normal  to  the 
edge. 

The  boundary  conditions  for  the  four  point  bend 
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specimen  are  (Figure  30b) : 

1.  °  ^(w.y)  -  0,  0  <  y  <  h, 

a  xx(v,y)  -  -Pi(y-h) , 

2.  0  xyCO.y)  -  0,  0  <  y  <  (h+d), 

o  xx(O.y)  “  Pi(y-h-d), 

3.  a  xy(x,0)  -  0,  0  <  X  <  w, 

0  yy(x,o  )  -  0,  0  <  X  <  a, 

v(x,0)  -  0,  a  <  X  <  w, 

4 •  C  yy (X , L)  •  C xyCx , L)  m  0 . 

In  the  finite  element  analysis,  the  boundary  condition 
c  xxCw.y^  “  -PS(y-ii)  was  replaced  by  u(w,h)  -  0  to  prevent 
rigid  body  displacement . 

For  the  three  point  bend  specimen,  the  boundary 
conditions  are  (Figure  30c): 

1.  "xy(w*y)  0  <  y  <  L, 

0  xx(w-y)  "  -PKy-0), 

2.  a  xy(0 ,y)  -  0,  0  <  y  <  h, 

axxC°-y)  "  PKy-h), 

3.  axy(x.°)  “0.  0  <  x  <  w, 

cyy(x*o)  “  o  <  x  «  a, 

v(x,0)  -  0,  a  <  x  <  v, 

4.  ayy(X,L)  “  0jfy(X,L)  “  0. 


43 


NADC-87 145-60 


The  boundary  condition  °xx(v,y)  -  -P5(y-0)  was  replaced  by 
u(v,0)  -  v(w,0)  *  0  in  the  finite  element  analysis  to 
prevent  rigid  body  displacement. 

The  nominal  bending  stress  for  the  three  and  four 
point  bend  specimens  is  calculated  from 

5 O  -  M0C/I  -  6M0/v2  (40) 

where  c-w/2,  I  is  the  area  moment  of  inertia  about  the  z- 
axis,  and  M0  -  p-h  for  the  three  point  and  M0  -  p-d  for 
the  four  point  specimen. 

This  problem  utilized  the  same  meshes  and  geometry  in 
the  applied  load  boundary  condition  comparison.  Figure  19. 
The  ASTM  [19]  dimensions  for  the  bending  specimens  are  a/w 
-  .5  and  h/w  -  2. 

The  normalized  stress  intensity  factor,  K]_/ “0(7ra)1/2 , 
is  shown  for  various  half-length  to  width  ratios  (h/w)  in 
Figure  31a, b  and  Table  10  for  the  bending  specimens.  The 
normalized  stress  Intensity  factors  for  the  four  point 
bend  and  pure  bend  specimen  are  approximately  the  same  for 
h/w  >  .8,  however  the  three  point  bend  specimen  is 
significantly  different  from  the  other  specimens  until 
h/w  2.  5.  The  normalized  stress  intensity  factor  of  1.518 
for  the  pure  bending  specimen  is  1.5  percent  different 
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from  the  theoretical  value  of  1.496  computed  from  [16], 
and  the  value  of  1.424  for  the  three  point  bend  specimen 
(for  h/w  -  2)  is  0.9  percent  different  from  the 
theoretical  value  of  1.411  computed  from  [16]. 

Another  study  was  performed  on  the  four  point  bend 
specimen  to  determine  the  effects  of  changing  the  distance 
between  the  point  load  and  the  roller  pin,  d  (Figure  29b). 
A  specimen  with  h/w  -  10,  a/w  -.5,  was  used.  The  position 
of  the  point  load  P  was  kept  constant  at  the  coordinates 
(0,9w)  while  the  position  of  the  pin  was  varied.  The 
results  in  Table  11  show  that  changing  d  had  no  effect  on 
the  normalized  stress  intensity  factor. 

Deformed  plots  of  the  bending  specimens  are  shown  in 
Figure  32a, b,c.  The  uniform  bend  and  three  point  bend 
specimens  show  more  rotation  than  the  four  point  bend 
specimen. 

Finally  the  normal  stress  distribution  in  front  of 
the  cracfc  tip,  ryy(x,0),  is  provided  in  Figures  33a, b,c 
and  Tables  12a, b,c  for  the  bend  specimens.  The  normal 
stress  distributions  for  the  pure  bend  specimen  and  the 

four  point  bend  specimen  are  very  similar,  while  the  three 

> 

point  bend  specimen  has  slightly  different  results.  As  h/w 
decreases,  the  normal  stress  distribution  for  the  pure 
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bend  and  four  point  specimens  changes  very  little,  however 
the  distribution  changes  significantly  for  the  three  point 
bend  specimen.  There  are  large  compressive  stresses  in  the 
three  point  bend  specimen  near  x/v  -  1.0  (near  the  pin) 
for  h/w  -  .5. 
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4.  CONCLUSION: 

Accurate  values  of  the  Mode  I  stress  intensity  factor 
were  computed  using  the  collapsed  quarter  point 
quadrilateral  elements,  however  the  mesh  around  the  crack 
tip  must  he  carefully  prepared  which  can  be  tedious.  These 
elements  used  with  the  direct  method  provide  a  relatively 
easy  way  to  compute  stress  intensity  factors  for  complex 
geometries  and  loading  conditions  where  handbook  solutions 
do  not  exist.  Also,  the  crack  opening  displacements  from 
the  finite  element  analysis  yield  more  accurate  stress 
intensity  values  than  the  stresses  in  front  of  the  crack 
tip. 


The  single  edge  cracked  specimen  under  the  flxeu 
grip,  displacement  controlled  loading  condition  had 
significantly  different  normalized  stress  intensity 
factors  than  the  point  and  distributed  loading  conditions 
for  half-length  to  width  ratios  (h/w)  less  than  twenty. 

The  normalized  stress  intensity  factors  for  the  point  and 
distributed  loading  conditions  were  exactly  tLe  same  after 
h/w  _>  1.25.  The  type  of  applied  loading  condition  effects 
the  stress  intensity  factor  around  the  crack  tip  for  this 
specimen. 
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The  normalized  stress  intensity  factor  for  the  three 
point  bend  specimen  was  significantly  different  from  the 
pure  bend  and  four  point  bend  specimen  for  h/w  <  5.  The 
normalized  stress  intensity  factors  for  the  pure  bend  and 
four  point  bend  specimens  were  approximately  the  same 
except  for  very  short  specimens.  The  four  point  bend 
configuration  produces  a  closer  approximation  to  a  pure 
bending  load  than  the  three  point  bend  configuration. 
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TABLES 
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Theta.  6  .  (degrees)  — °  6J,Z  ± 


0.00 

3.99 

3.75 

3.95 

7.50 

3.86 

11.25 

3.69 

15.00 

3.49 

18.75 

3.22 

22.50 

2.92 

26.25 

2.59 

30.00 

2.22 

33.75 

1.86 

37.50 

1.47 

41.25 

1.11 

45.00 

.77 

48.75 

.44 

52.50 

.  12 

56.25 

-.20 

60.00 

-.46 

63.75 

-.70 

67.50 

-.89 

71.25 

-1.05 

75.00 

-1.16 

78.75 

-1.25 

82.50 

-1.30 

86.25 

-1.33 

90.00 

-1.34 

Table  1.  Tangential  stress  ratio,  " -  z  t  at  nodal 

points  about  the  boundary  of  the  large  central 
hole . 
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Thet a .0. (degrees j 


o 


-90 

-.11 

-81 

-.11 

-72 

-.11 

-63 

-.05 

-54 

.10 

-45 

.28 

-36 

.55 

-27 

.86 

-18 

1.18 

-9 

1.51 

0 

1.76 

9 

1.84 

18 

1.79 

27 

1.60 

36 

1.26 

45 

.97 

54 

.63 

63 

.  17 

72 

-.15 

81 

-.38 

90 

-.44 

Table  2.  Tangential  stress  ratio,  "  -/  0  ,  at  nodal 

points  about  the  boundary  of  the  small  holes. 
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Crack 


Length, a( in) 

*l/°o 

Ki/^oOra) 

.01 

.  704 

3.973 

.02 

.949 

3.784 

.04 

1.226 

3.458 

,06 

1.405 

3.235 

.08 

1.546 

3.083 

.10 

1.671 

2.981 

.125 

1.823 

2.910 

.15 

1.979 

2.883 

.20 

2.342 

2.954 

.25 

2.846 

3.211 

.30 

3.703 

3.814 

.325 

4.510 

4.464 

.35 

6.298 

6.006 

1/2 


Table  3.  Normalized  stress  intensity  factor  vs.  crack 
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Kl/a0(7Ta)1/2 


h/w 

Point  Load 

Distributed  Load 

Fixed  ( 

.  1 

5.137 

14.871 

.455 

.2 

3.635 

6.254 

.645 

.3 

3.013 

4.097 

.778 

.4 

2.775 

3.344 

.873 

.5 

2.720 

3.035 

.959 

.6 

2.733 

2.914 

1.037 

.8 

2.788 

2.646 

1.171 

1.0 

2.823 

2.836 

1.287 

1.25 

2.836 

2.836 

1.415 

1.5 

2.839 

2.836 

1.524 

3 

- 

- 

1.940 

5 

- 

- 

2.208 

10 

- 

- 

2.470 

15 

- 

- 

2.581 

20 

- 

- 

2.637 

30 

- 

- 

2.702 

40 

- 

- 

2.733 

100 

2.836 

2.836 

2.793 

Table  4.  Normalized  stress  intensity  factors  vs.  half- 
length,  to  width  ratio  (h/w)  for  the  single 
edge  crached  specimen  under  various  applied 
loading  conditions. 
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AU£  /(  a  o  h/E) 


h/w 

Point .Load 

Distributed  Load 

Fired  Grip 

.  1 

10.5862 

54.9415 

.8345 

.2 

5 . 3400 

15.1480 

.7984 

.3 

4.1126 

7 . 9745 

.7555 

.4 

3.9190 

5 . 8678 

.7318 

.5 

4.0403 

5.1455 

.7414 

.6 

4.2541 

4.9094 

.7701 

.6 

4.6752 

4.8923 

.8554 

1.0 

4.9536 

5.0019 

.9548 

1.25 

5.1296 

5.1191 

1.0780 

1.5 

5.2134 

5.2001 

1.1915 

3 

- 

- 

1.6671 

5 

- 

- 

1 . 9968 

10 

- 

- 

2.3332 

15 

- 

- 

2.4734 

20 

- 

- 

2 . 5468 

30 

- 

- 

2.6283 

40 

- 

- 

2.6695 

100 

5 . 5984 

5.5984 

2.7477 

Table  5.  Normalized  transverse  relative  displacement 
(J^u^/c  0  h/E  )  vs  half-length  to  width, 
ratio  (h/v)  for  the  single  edge  cracked 
test  specimen  under  various  loading 
conditions . 
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Point  Load 
°  yy^* 1 0) /  0  o 


x/v 

h/w  -  100 

h/w  -  1.5 

.5 

54.133 

54.184 

.5125 

12.116 

12.124 

.55 

8.842 

8.847 

.575 

4.372 

4.371 

.600 

3.506 

3.505 

.625 

2.920 

2.917 

.65 

2.423 

2.420 

.675 

2.012 

2.009 

.700 

1.639 

1.635 

.725 

1.313 

1.309 

.750 

.  977 

.974 

.8125 

.229 

.227 

.875 

-.581 

-.581 

.9375 

-1.500 

-1.497 

.0 

-2.731 

-2.726 

Table  6a.  Normal  stress  in  front  of  the  crack  tip. 

c  jryCx.O).  for  the  single  edge  cracked 
specimen  under  a  point  load. 
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Point  Load 

0  yy  • o) /  c  q 


x/w 


h/W  -  .  1 


.5 

.50623 

.525 

.5375 

.550 

.5625 

.5750 

.5875 

.600 

.625 

.650 

.675 

.700 

.725 

.750 

.775 

.800 

.825 

.850 

.875 

.900 

.925 

.950 

.975 

1.0 


138.294 

31.251 

22.511 

9.290 

5.841 

3.365 

1.593 

.280 

-.504 

-1.230 

-1.252 

-.970 

-.690 

-.421 

-.257 

-.131 

-.067 

-.026 

-.008 

.001 

.003 

.004 

.003 

.003 

.004 


TaJbie  6a  (continued). 
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Distributed  Load 


a  yy(x •  ^  0  o 


x/v 

h/w  -  100 

b/v  -  1 

.5 

54.133 

54.134 

.5125 

12.116 

12.116 

.55 

8.842 

8.842 

.575 

4.372 

4.372 

.600 

3.506 

3.506 

.625 

2.920 

2.920 

.65 

2.423 

2.423 

.675 

2.012 

2.012 

.700 

1.639 

1.638 

.725 

1.313 

1.313 

.750 

.977 

.977 

.8125 

.229 

.229 

.875 

-.581 

-.581 

.9375 

-1.50 

-1.500 

1.0 

-2.731 

-2.731 

Table  6b.  Normal  stress  in  front  of  the  crack  tip, 

-  yy (x,0),  for  tie  single  edge  cracked 
specimen  under  a  uniformly  distributed  load. 
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Distributed  Load 


0  yy(x,0)/  0  o 


x/v 

b/w  -  . 1 

.5 

426 . 646 

.50625 

66.295 

.525 

38.586 

.5375 

.705 

.550 

-5.909 

.5625 

-8.981 

.5750 

-10.511 

.5875 

-10.575 

.600 

-10.104 

.625 

-7.606 

.650 

-5.155 

.675 

-2 . 772 

.700 

-1.334 

.725 

-.208 

.750 

.374 

.775 

.751 

.800 

.913 

.825 

.999 

.650 

1.025 

.675 

1.030 

.900 

1.025 

.925 

1.018 

.950 

1.013 

.975 

1.009 

1.0 

1.010 

Table  6b  (continued). 
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Fixed  Grip 


0  yy (X , 0 ) / 0  q 


x/w 

h/v  -  100 

h/v  -  1.5 

.5 

53.324 

28 . 778 

.5125 

11.945 

6.770 

.55 

8.722 

5.066 

.575 

4.322 

2.807 

.600 

3.472 

2.419 

.625 

2.896 

2.175 

.65 

2.409 

1.975 

.675 

2.006 

1.821 

.700 

1.640 

1.683 

.725 

1.321 

1.571 

.750 

.992 

1.458 

.8125 

.260 

1.213 

.875 

-.532 

.945 

.9375 

-1.432 

.625 

1.0 

-2.639 

.160 

Table  6c.  Normal  stress  in  front  of  the  crack  tip, 
-yy(x,0),  for  the  single  edge  cracked 
specimen  under  fixed  grip-displacement 
controlled  conditions. 
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Fixed  Grip 
o  yy(x > 0)  f  a  o 


x/w 


h./v  -  .  1 


.5 

.50625 
.525 
.5375 
.  550 
.5625 
.5750 
.5875 
.600 
.625 
.650 
.675 
.700 
.725 
.750 
.775 
.800 
.825 
.850 
.875 
.900 
.925 
.950 
.975 
1.0 


11.566 

3.309 

2.684 

1.973 

1.923 

1.920 

1.920 

1.932 

1.941 

1.956 

1.966 

1.970 
1.972 

1.971 

1.972 
1.969 
1.969 
1.965 
1.963 
1.958 
1.952 
1.913 
1.847 
1.708 
1.408 


Tafcle  6c  (continued). 
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Singled  Edged 
Cracked  Specimen 
(Fixed  Grip) 

*1 _ 

Ev0/[h(l-y2)]1/2 


Infinitely  Vide 
Strip  With  A 
Fixed  End 


Ev0/  [hCl-i'2)  ] 1/2 


.  1 
.2 
.3 
.4 
.5 
.6 
.8 
1.0 
1 . 25 
1.5 
3 
5 

10 

15 

20 

30 

40 

100 


.  9653 
.978? 
.9802 
.9699 
.9689 
.  9688 
.9742 
.9834 
.9952 
1 . 0043 
1.0048 
.9535 
.8200 
.7234 
.6531 
.  5571 
.4935 
.3257 


Table  7.  A  comparison  between  the  single  edge  cracked 
test  specimen  under  fixed  grip  conditions  and 
tie  infinitely  wide  strip  solution  with  a 
fixed  end. 
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h/V  m  .  1 


:/w 

A 

B 

x/w 

A 

B 

). 

.01 

o 

o 

.6 

1.94 

-.05 

.025 

.00 

.625 

1.95 

-.03 

.050 

.650 

1.95 

.00 

.075 

.675 

1.95 

.01 

.100 

.700 

1.95 

.02 

.125 

1 

t 

.725 

1.94 

.03 

.150 

.00 

.750 

1.93 

.05 

.175 

-.01 

.775 

1.92 

.07 

.200 

.800 

1.91 

.09 

.225 

.825 

1.89 

.  11 

.250 

> 

f 

.850 

1.87 

.14 

.275 

-.01 

.875 

1.85 

.  18 

.300 

-.02 

.900 

1.84 

.23 

.325 

\ 

f 

-.03 

.925 

1.75 

.27 

.350 

.00 

-.05 

.950 

1.72 

.33 

.375 

.03 

-.09 

.975 

2.05 

.  53 

.400 

.07 

-.15 

1.0 

2.85 

.79 

.425  .20  -.25 
.450  .44  -.34 
.475  .84  -.39 
.500  1.27  -.40 
.525  1.61  -.30 
.550  1.84  -.19 
.575  1.92  -.12 


A  “  yy(X,ll)/rQ 
B  “  xy  C  x ,  h )  /  r  g 

Table  8.  Normal  stress  C-'yy)  and  shear  stress  (  ;xy^ 
along  the  fixed  grip  end  of  the  single  edge 
cracked  speciaen. 
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h./V  -  .5 


h/v  -  1.5 


X/v  A  B 


A  B 


). 

.01 

.0625 

-.01 

.125 

.04 

.  1875 

.17 

.2500 

.33 

.3125 

.49 

.3750 

.68 

.4375 

.88 

.5000 

1.07 

.5625 

1.24 

.6250 

1.36 

.6875 

1.48 

.7500 

1.57 

.8125 

1.58 

.8750 

1.66 

.9375 

2.12 

.0 

2.95 

-.03 
-.11 
-.15 
-.17 
-.20 
-.21 
-.22 
-.19 
-.15 
-.09 
-.02 
.06 
.  15 
.21 
.26 
.50 
.77 


.44 

-.18 

.39 

-.17 

.37 

-.15 

.45 

-.  14 

.54 

-.13 

.63 

-.12 

.73 

-.10 

.83 

-.08 

.93 

-.06 

1.03 

-.03 

1.14 

.00 

1.25 

.04 

1,36 

.08 

1.48 

.14 

1.62 

.22 

2.02 

.37 

2.45 

.52 

A  "  '  yy(X,ll)/-Q 
B  -  -xy(X,h)/;0 

Table  8 .  (coat . ) 
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h/V  -  10.  h/V  -  100. 


x/v 

A 

B 

A 

B 

1. 

1.022 

-.244 

1.26 

-.27 

.0625 

.939 

-.206 

1.14 

-.22 

.  1250 

.859 

-.166 

1.02 

-.17 

.  1875 

.859 

-.138 

.99 

-.14 

.2500 

.863 

-.110 

.97 

-.  11 

.3125 

.881 

-.085 

.96 

-.08 

.3750 

.  901 

-.061 

.95 

-.05 

.4375 

.927 

-.037 

.95 

-.03 

.5000 

.954 

-.013 

.95 

.00 

.5625 

.985 

.012 

.96 

.02 

.6250 

1.018 

.038 

.97 

.05 

.6875 

1.057 

.069 

.98 

.08 

.7500 

1.099 

.  101 

.99 

.10 

.8125 

1.159 

.  142 

1.03 

.  14 

.8750 

1.223 

.184 

1.06 

.18 

.9375 

1.386 

.244 

1.19 

.23 

.0 

1.552 

.305 

1.32 

.28 

A  “  -  yy  C  X  ,  h  )  /  ”  q 
B  -  c  0 

Table  8 .  (cont . ) 
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r/w 

0. 

.0625 

.1250 

.1675 

.2500 

.3125 

.3750 

.4375 

.5000 


SIDE  A 
_ A _ 

.2154 

.1956 

.1705 

.1400 

.1093 

.0017 

.0536 

.0276 

.0017 


SIDE  B 

_ B_ 

.2107 
.1900 
.1656 
.1364 
.  1053 
.0701 
.0502 
.0243 
.0017 


A  -  _cXy(r.h)/^yy(r.lx)  For  Side  A 
B  -  axy(r  ,h)/3yy(r.h)  For  Side  B 


Table  9.  Shear  stress  to  normal  stress  ratio  at  the 
fixed  grip  end  of  the  single  edge  cracked 
specimen  (h/v  -  100). 
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Kj./  a  0(Tr*')l/2 


Pure 

4  Point 

3  Point 

h/w 

Beading 

Bending 

Sending 

.5 

1.650 

1.555 

1.158 

.6 

1.570 

1.546 

1.217 

.8 

1.525 

1.530 

1.289 

1.0 

1.518 

1.522 

1.332 

1.8 

1.518 

1.392 

2.0 

1.517 

1.424 

5.0 

M 

f 

1.517 

1.480 

10.0 

1.518 

1.517 

1.499 

Crack  length  to  width  ratio,  aA  -  .5 


Table  10.  A  comparison  of  normalized  stress  intensity 

factors  for  single  edge  cracked  bend  specimens 
for  various  half-length  to  width  (h/w)  ratios. 
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d/W 


Single  Edge  Cracked 
4  Point  Bead  Specimen 

Kl/a  0(7Ta)1/2 


.25 

.5 

1.0 

1.5 

2.0 

7.0 


1.318 

1.517 

1.5 18 

1 

1.518 


Crack  length  to  width  ratio,  a/w  -  .5 
Specimen  length  to  width  ratio,  h/w  -  io 


Table  ii. 


®tr!ss  inteasity  factor  of  the 
four  point  bend  specimen  for  various 

distances  between  the  point  and  the  roller 
pin,  d,  (see  Figure  29b). 
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X/V 


.5 

.  5125 
.55 
.575 
.6 

.625 

.65 

.675 

.7 

.725 

.75 

.8125 

.875 

.9375 

1.0 


Pure  Bend  Specimen 

0  yy^*’ o 
h/v-10  h/v-1.5 


29.345 

6.188 

4.371 

1.812 

1.259 

.863 

.518 

.222 

-.052 

-.300 

-.557 

-1.140 

-1.766 

-2.460 

-3.345 


29.345 
6. 188 
4.371 
1.812 
1.259 
.863 
.518 
.222 
-.052 
-.300 
-.557 
-1.140 
-1.766 
-2.460 
-3.345 


h/w  - 


32.103 

6.558 

4.558 
1.735 
1.129 

.707 
.351 
.056 
-.206 
-.438 
-.669 
-1.176 
-1 . 726 
-2.365 
-3.277 


Ta*le  12a.  Normal  stress  in  front  of  the  crack  tip 

-yy(x.O),  for  the  single  edge  cracked 
specimen  under  pure  binding8  Crac3ced 
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4  Point  Bend  Specimen 


OyyCx.0 )/o0 


x/v 

h/v-10 

h/V-1.5 

h/v-  . 5 

.5 

29.345 

29.346 

30.109 

.5125 

6.188 

6.187 

6.311 

.55 

4.371 

4.371 

4.444 

.575 

1.812 

1.812 

1.808 

.6 

1.259 

1.259 

1.236 

.625 

.863 

.862 

.827 

.65 

.518 

.518 

.472 

.675 

.222 

.222 

.168 

.7 

-.052 

-.052 

-.111 

.725 

-.300 

-.300 

-.359 

.75 

-.557 

-.557 

-.614 

.8125 

-1.140 

-1.140 

-1.180 

.875 

-1.766 

-1.765 

-1.767 

.9375 

-2.460 

-2.460 

-2.414 

1.0 

-3.345 

-3 . 346 

-3.257 

Table  12b.  Normal  stress  in  front  of  the  crack  tip, 
Cyy(x.O),  for  the  single  edge  cracked 
specimen  under  four  point  bending. 
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I  -TW  wv 


Figure  2a.  A  2D  solid  with  a  flat  crack, 


Y 


Figure  2b.  Stress  state  around  the  crack  tip. 


73 


NADC-871 45-60 


Y 


Figure  6a.  Boundary  conditions  for  the  centrally  cracked 
plate. 


77 


NADC-87 145-60 


Figure  10 
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Finite  element  mesh  of  fatigue  specimen 
without  a  crack. 
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Figure  14a.  Finite  element  mesh  of  the  fatigue  specimen 
with  a  crack. 
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Region  A 


Figure  14b.  Enlarged  view  of  lower  mesh  region. 
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Normalized  stress  intensity  factor  vs.  crack  length  for 
fatigue  test  specimen. 
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Stress  intensity  factor  normalised  with  respect  to  OJ,  vs. 
track  length  for  fatigue  test  specimen. 
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K -  a  - * 


Figure  19a.  Finite  element  mesh  for  the  single  edge 
cracked  specimen,  h/w  =  .5. 
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r -  a  - H 


Figure  19b. 


Finite  element  mesh  for  the  single  edge 
cracked  specimen,  h/w  =  1.5. 
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Figure  19c.  Finite  element  mesh  for  the  single  edge 
cracked  specimen,  h/w  *  20. 
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loading  conditions. 


FINITE  ELEMENT  RESULTS 
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Figure  21a.  Deformed  plot  of  the  point  loaded  single  edge 
cracked  specimen,  h/w  *  1.5. 
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Figure  21b.  Deformed  plot  of  the  uniformly  distributed  loaded 
single  edge  cracked  specimen,  h/w  -  1.5. 
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Figure  21c.  Deformed  plot  of  the  single  edge  cracked  specimen 
under  displacement  controlled  -  fixed  grip 
conditions.  aq 
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Figure  21d.  Deformed  plot  of  the  point  loaded  single  edge 
cracked  specimen,  h/w  *  .1. 
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Figure  2le.  Oeformed  plot  of  the  uniformly  distributed  loaded 
single  edge  cracked  specimen,  h/w  *  .1. 
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Figure  21f.  Deformed  plot  of  the  single  edge  cracked  specimen 
under  displacement  controlled  -  fixed  grip 
conditions. 


102 


Figure  22a.  Normalized  transverse  relative  displacement  (  — - 1 —  )  vs  half-length 

to  width  ration  (h/w)  for  the  single  edge  cracked  specimen  under 
various  applied  loading  conditions. 


Figure  22b.  Norma  1 1  zed  transverse  relative  displacement  (  ^ut  )  vs  half-length 
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cracked  specimen  under  a  uniformly  distributed  load 


edge  cracked  specimen  under  fixed  grip  -  displacement  controlled 
conditions. 


A  comparison  betweer  the  single  edge  cracked  specimen  under 
fixed  grip  conditions  and  the  infinitely  wide  strip  with  a 
fixed  end. 
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Figure  26b.  A  comparison  between  the  single  edge  cracked  specimen  under  fixed 
grip  conditions  and  the  infinitely  wide  strip  with  a  fixed  end. 
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of  the  single  edge  cracked  specimen,  h/w 


NADC-871 45-60 


114 


Figure  27d.  Normal  stress  (0^y)  and  shear  stress  (CTXy)  along  the  fixed  grip 
end  of  the  single  edge  cracked  specimen,  h/w  =  10. 
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Figure  27>e.  Normal  stress  (CTyy)  and  shear  stress  (O^y)  along  the  fixed  grip 
end  of  the  single  edge  cracked  specimen,  h/w  =  100. 
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Figure  29b.  Four  point  bend  specimen- 
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Y 
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Figure  29c.  Three  point  bend  specimen. 
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(a)  uniform  bend  specimen,  (b)  four  point  bend  specimen,  and  (c)  three  point 
bend  specimen. 


Figure  31a.  Normalized  stress  intensity  factor  vs.  half-length  to  width  ratio 
(h/w)  for  single  edge  cracked  bend  specimens. 
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Normalized  stress  intensity  factor  vs.  half-length  to  width  ratio 
(h/w)  for  single  edge  cracked  bend  specimens. 


specimen 


--•r - 4 . .  ■ 

Figure  32b.  Deformed  plot  of  a  single  edge  cracked  specimen 
under  four  point  bending. 
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•l*. 


specimen 


Figure  33a.  formal  stress  in  front  of  the  crack  tip,  0^y(x,o),  for  the  single 
edge  cracked  specimen  under  pure  bending. 
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O 


Normal  stress  in  front  of  the  crack  tip,  O^y(x,o),  for  the  single 
edge  cracked  specimen  under  four  point  bending. 


1  stress  in  front  of  the  crack  tip,  (fyy(x,o),  for  the  single 
cracked  specimen  under  three  point  bending. 
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